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NEAR-IR POLARIMETRY AROUND 30 DORADUS: 
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ABSTRACT 

A 20' x 20' region around 30 Doradus in the Large Magellanic Cloud (LMC) is observed and ana- 
lyzed in the near-infrared. We obtain polarimetry data in the J, H, and K s bands using the SIRIUS 
polarimeter SIRPOL at the Infrared Survey Facility 1.4 m telescope. We measure the Stokes parame- 
ters of 2562 point-like sources to derive the degree of polarization and the polarization position angles. 
We discuss the statistics of the groups classified by color-magnitude diagram and proper motions of 
the sources, in order to separate the Galactic foreground sources from those present in the LMC. We 
notice that groups classified by the proper motion data show a tendency towards different polarimetric 
properties. 
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1. INTRODUCTION 

The role of magnetic fields is relevant in star for- 
mation processes and in galaxy evolution science (Ke- 
pley 2009; Schmidt 1970, 1976; Mathewson and Ford 
1970; Wielebinski 1995). Measuring the orientation of 
dust grains is a way to investigate the physical prop- 
erties of the magnetic field, since the magnetic field 
makes dust grains to align, due to rapid precession en- 
forced from the substantial magnetic moment (Davis & 
Greenstein 1951; Dolginov & Mytrophanov 1976). In 
star-forming molecular clouds, dust grains tend to be 
aligned with their long axes perpendicular to the mag- 
netic field (Lazarian 2006); dichroic polarization in the 
optical and near-IR wavelengths can measure the di- 
rection of those fields (Davis & Greenstein 1951). In 
addition, the statistical pattern of the polarization po- 
sition angles allows one to quantify the magnetic field 
strength (Ostriker et al. 2001; Kwon et al. 2010). 
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The Large Magellanic Cloud (LMC) is one of the 
nearest external galaxies. The interstellar medium of 
the LMC shows quite different physical and chemical 
properties with respect to that of our Galaxy, e.g., 
lower metallicity and different dust properties (Pak et 
al. 1998; Nakajima et al. 2007; Kim et al. 2010). The 
LMC has revealed large-scale magnetic field structures 
(Gaensler et al. 2005) from optical polarization ob- 
servations (Schmidt 1970, 1976; Mathewson and Ford 
1970; Wayte 1990). These features indicate that, in 
the west side of the LMC, the magnetic field is point- 
ing towards the Small Magellanic Cloud. On the other 
hand, the magnetic fields around 30 Doradus appear to 
be complex structures. In fact, 30 Doradus is a unique, 
giant star formation complex in the LMC and contains 
massive clusters (Walborn & Blades 1997; Brandner et 
al. 200 l;Maercker& Burton 2005 ;Scowenetal. 2009). 
Multi-frequency radio continuum surveys (Haynes et al. 
1991) show the presence of filamentary magnetic fields 
in the LMC. 

For this study, we observed near-infrared (NIR) 
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sources in star forming regions (~39'x69') in the LMC. 
The number of detected sources, however, is quite 
large. In this paper, we only present preliminary results 
of selected regions around 30 Doradus (~20'x20'). Sec- 
tion 2 describes the observations and the data reduc- 
tion processes, including photometry and the deriva- 
tion of the Stokes parameters. Section 3 includes the 
comparison of our data with previous observations by 
Nakajima et al (2007). In Section 4, we discuss the 
methods to extract Galactic foreground sources using 
the color-magnitude diagram and proper motions. 

2. OBSERVATIONS AND DATA REDUC- 
TION 

The observations were carried out at the Infrared 
Survey Facility (IRSF) 1.4 m telescope in the South 
African Astronomical Obsevatory on December 25 and 
30, 2008. Using the infrared camera SIRIUS (Na- 
gayama et al. 2003) and the polarimcter SIRPOL 
(Kandori et al. 2006), we observed regions around 30 
Doradus in J (1.25 /zm), H (1.63 /zm), and K s (2.14 
/an) bands. This system enables us to obtain a wide- 
field (7.'7 x 7.'7) polarimetric image with a scale of Of! 45 
pixel . 

We performed four 20-second exposures at wave- 
plate angles of 0°, 45°, 22°. 5, 67°. 5, and at 10 dithered 
positions for each field position. The total integration 
time was 200 seconds per wave-plate angle. The typical 
seeing size was ~l'/3 in the J band, but we experienced 
unstable seeing conditions during a night. Among the 
observed 5x9 fields (^39'x69') we analyzed the central 
3x3 fields (~20'x20') and illustrate here our initial in- 
vestigations (see Fig. 1). Table 1 lists the observation 
log of the data presented in this paper. 

We used the IRAF* (Image Reduction & Analy- 
sis Facility) software package to reduce the data, e.g, 
dark-field subtraction, flat-field correction and median 
sky subtraction. In addition, the data for each wave- 
plate position (Iq, J45, I22.5 and Ie7.s) were processed 
as described in Kandori et al. (2006). For source de- 
tection and photometry on the Stokes I images, we 
used the IRAF daophot package (Stetson 1987). We 
set the aperture radius to 8 pixels and the sky annulus 
to 13 pixels, and detected point sources having peak 
intensity greater than 10a above the local sky back- 
ground. The pixel coordinates of the detected sources 
were matched to the celestial coordinates of their coun- 
terparts in the 2MASS^ Point Source Catalogue. The 
magnitude and color of our photometry were trans- 
formed into the 2MASS system using: 



"IRAF is distributed by the US National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 

+ The Two Micron All Sky Survey(2MASS) is a joint project of 
the University of Massachusetts and the Infrared Processing and 
Analysis Center / California Institute of Technology. 



m 2 MASS = m IRSF + o-i X C IRSF + (1) 
C2MASS — Ct2 X ClRSF + 02, (2) 

where ttiirsf is the instrumental magnitude from the 
IRSF images and rri2MASs is the magnitude from the 
2MASS All Sky Point Source Catalogue. 

We calculated the Stokes parameters I, Q, and U, 
the polarization degree P, and the polarization posi- 
tion angle 0, as follows: 

/ = (Jo + ^22.5 + h 5 + / 6 7.5)/2, (3) 

Q = I -h 5 , (4) 

U = 722.5 - ^67.5, (5) 

P. . v^E, (6) 

= ^ arctan ^. (7) 

Since the initial polarization degree P is a positive 
quantity, the derived P Q values tend to be overesti- 
mated, especially for low signal-to- noise ratio (S/N) 
sources. To correct for the bias, we use the following 
equation: 

p = VpF~sp!, (8) 

where 5P is the error in P (Wardle & Kronberg 1974). 
The polarization efficiencies of SIRPOL (95.5, 96.3, and 
98.5% at J, H, and K s bands, respectively) were also 
used for correcting the polarization degree (Kandori et 
al. 2006). 

3. RESULTS 

From the data reduction process, we first derive po- 
larization degrees and angles of 2562 sources, and ex- 
tract only those with P > 3cr p. (op indicates the error 
of a polarization value.) The numbers of the selected 
stars are 736, 784, and 427 in J, H, and K s bands, 
respectively. Fig. 2 shows a polarization vector map in 
the reduced 9 fields. The statistical distribution of the 
polarization degrees and angles is shown in Fig. 3. 

To confirm our results, we compare our data with 
previous observations. Nakajima et al. (2007), here- 
after N07, had observed the central 7'.7 x 7.'7 (one field) 
regions of 30 Doradus using the same SIRIUS/SIRPOL 
system. Our photometry in the 30 Doradus is limited 
by the nebula features, while N07 detected many point 
sources with long exposures, i.e., a total of 1480 seconds 
per wave-plate angle. 

The numbers of the matched sources in both data 
are 33, 36, and 7, at J, H, and K s bands, respectively. 
With these sources, we compared the polarization de- 
gree and angle. Figure 4 shows that our values of the 
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Table 1. 

Observation Log 



Datc(LT) 


Field name 


C*oJ2000.0 


OoJ2000.0 


Seeing in J band ["] 


2008 Dec 25 


n4 


05 37 19.86 


-69 06 03.3 


1.4 




n5 


05 38 34.48 


-69 06 02.2 


1.4 




n6 


05 39 48.78 


-69 06 02.6 


1.3 




n7 


05 39 48.17 


-69 12 38.0 


1.3 




n8 


05 38 32.91 


-69 12 34.0 


1.2 




n9 


05 37 17.79 


-69 12 32.4 


1.3 


2008 Dec 30 


nl 


05 39 47.92 


-68 59 09.9 


1.2 




n2 


05 38 33.53 


-68 59 07.7 


1.3 




n3 


05 37 19.27 


-68 59 03.9 


1.3 



Table 2. 

Statistical Properties of the Polarization angle 





J a 




K s a 


All bands" 


All sources 


79 ± 47 


82 ± 48 


85 ± 46 


74 ± 24 


GroupB 
GroupEFG 


75 ± 49 
81 ± 47 


83 ± 53 
85 ± 47 


75 ± 43 
87 ± 49 


76 ± 11 
73 ± 21 


[l > 2(7^ 
fJ. < 2cr M 


90 ± 45 
69 ± 46 


102 ± 51 

82 ± 48 


97 ± 44 
83 ± 46 


96 ± 19 
65 ± 24 



Note. — a Polarization angle in degree units. The central 
angle and the sigma values are from obtained Gaussian fits to 
the distribution functions. 
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polarization angles and those from N07 are somehow 
consistent. However, the polarization degree values 
in our works are 50% larger than those in N07. We 
suspect that this discrepancy is caused by the under- 
estimated SP values in our low S/N data. More ob- 
servations and analysis is required to confirm the cali- 
bration of the polarization degree values. 

4. DISCUSSION 

Source classification is important to separate Galac- 
tic foreground sources from those in the LMC, because 
they can contaminate polarimetry by their intrinsic po- 
larization. In this section, we applied two methods to 
test the separation reliabilities. 

4.1 Color-Magnitude Diagram 

The typical classification method is to use color- 
magnitude diagrams (hereafter CM diagrams). Niko- 
laev & Weinberg (2000), hereafter NW00, classified the 
LMC stars using 2MASS data as follows: 

• Group A: blue supergiants and O dwarfs in the LMC; 

• Group B: Galactic disk dwarfs with spectral classes 
of F— K and Galactic giants with spectral classes of 
F-G; 

• Group D: Galactic disk dwarfs with spectral classes 
of G— M, red-giant branch (RGB), and early asymp- 
totic giant branch (E-AGB) stars in the LMC; 

• Group EFG: tip of the RGB, E-AGB, and O-rich 
AGB stars in the LMC ; 

• Group I: super giants and the red clump in the LMC; 

• Group JK: C-rich AGB and dusty AGB stars in the 
LMC. 

The 2MASS sources, however, are limited in itik s 
< 14 mag. On the other hand, our photometry data 
cover up to rriK a < 16 mag and most detected sources 
are populated under Group D in NW00. Kato et al. 
(2007) used the criteria of NW00, but did not clas- 
sify the sources which distributed under Group D. We 
followed the same criteria as Kato et al. (2007) to sep- 
arate sources. Figure 5 shows our five groups. 

Galactic foreground stars are mostly distributed 
above the brighter region {rriK s < 11) of Group A 
according to NW00. However, our data do not con- 
tain the detected sources in that region. The detected 
foreground sources are distributed in Group B; the 
Galactic giants with spectral classes of F— G have a 
small population (< 5%) within this group. Most 
RGB and AGB stars in the LMC arc in Group EFG. 
Note that Group A and Group I can be confused 
with the Galactic foreground sources (Group B) be- 
cause of their uncertain boundary lines. Group JK 
have dusty AGB stars which may have their own po- 
larization by the shrouded dust shells (Parthasarathy 
& Jain 1993). Therefore, Group A, Group I and Group 
JK are excluded from the samples when comparing the 
foreground and the main LMC sources. 



We classify Group B as foreground sources and 
Group EFG as the main LMC stars. We extrapolate 
the boundaries of the criteria regions up to m Ks < 16 
mag and plot histograms for the sources in Group B 
and Group EFG. Figure 6 and 7 show the distribu- 
tions of polarization degrees and angles for both Group 
B and Group EFG. In Figure 6, the polarization de- 
gree distribution of Group EFG shows a large pop- 
ulation of high polarization values (P > 10%) which 
comes from the intrinsic polarization of the RGB and 
AGB stars. On the other hand, the polarization angle 
distributions show similar tendency. See the statistical 
distributions in Figure 7 and Table 2. Because the cri- 
teria of Group B and Group EFG are overlapped in 
14 < m Ks < 16, we expect that many sources may not 
be classified properly. 

4.2 Proper Motion 

Vieira et al. (2010) investigated the proper motion 
of the Magellanic Clouds from the resulting catalog 
of 1.4 million objects in the Southern Proper Motion 
(SPM) program. They used the proper motion of 3822 
and 964 sources near the clouds to derive the mean 
absolute proper motions of the LMC and of the SMC. 

We used this catalog and examined proper motions 
to separate the foreground sources which would be con- 
tained in our observed regions in the LMC. A total 
of 805 sources were matched with SPM sources in the 
JHK S bands. Vieira et al. (2010) discarded areas of 
high stellar density, at 4 h 40 m 58.8 s < a < 5 /l 40 m 12 s 
and —68° < 8 < —71°, to avoid the misidentified mea- 
surement. Since our observed regions are located on 
the edge of the discarded area, some of our data may 
not be reliable. 

In order to reduce the uncertainty in the data, we 
used data which had proper motions with u > 2cr M to 
select Galactic foreground stars. Fig. 8 shows the dis- 
tribution of the sources according to the S/N of the 
proper motion. Based on the existence of proper mo- 
tion, we examined the difference between the LMC 
sources and the Galactic foreground ones. 

Figure 9 shows that the histograms of the polariza- 
tion degrees are similar between the two groups, while 
those of the CM diagram classifications in Figure 6 are 
very different. We expect that the difference in Figure 
6 is partly from the intrinsic differences of the spec- 
tral types. As evident from Figure 10, the polarization 
angle distributions are slightly different. The Galactic 
foreground group with proper motions shows almost 
equally distributed angles with a slight peak at 93°, 
while that of the LMC group shows a peak at 75°. 

We also examined the possibility that the Galactic 
foreground stars may show some polarization proper- 
ties different from the main LMC stars by their intrin- 
sic polarization or the effect for the Galactic magnetic 
fields. Previous observations (Mathewson & Ford 1970; 
Isserstedt & Rcinhardt 1976; Schmidt 1976) indicated 
that the Galactic magnetic fields at the line of sight of 
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the LMC can influence the polarization properties of 
the foreground stars. According to those papers, the 
polarization properties of the Galactic sources in front 
of the LMC have a 0.3% polarization degree and 30°of 
polarization angle. However, these results may not be 
applied to ours because our detection limits are too low 
to measure the dichroic polarization of the Galaxy. 

5. SUMMARY 

We conducted imaging polarimetry of the 3x3 fields 
(~20'x20') around 30 Doradus. Because of unstable 
seeing condition, our aperture radius for photometry 
and polarimetry has been set 2-3 times larger than the 
usual aperture radius. We used some previous data 
to confirm our results. Comparing with the sources 
of N07, we found no significant differences in our po- 
larimetry data. 

We applied several methods to remove the Galac- 
tic foreground sources which can contaminate the po- 
larimetry result by their intrinsic polarization. Using 
the CM diagram, we selected the group which contains 
the foreground sources. The main LMC sources show 
larger polarization degrees than the foreground sources. 

We also selected the Galactic foreground sources us- 
ing proper motion data. The foreground sources show 
high values of proper motion, contrary to the main 
LMC sources. 

Comparing different proper motions, we can sepa- 
rate out the Galactic foreground sources from the main 
LMC sources. They showed different distributions for 
polarization degrees and angles with respect to the 
main LMC sources. 
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Fig. 1. — Observed area of the LMC on an optical image offered by Mr. Kamiya. The yellow box shows the total 5x9 
survey field by SIRPOL on 2008 December 25 and 30 and the red box is the 3x3 fields around 30 Doradus. 
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Fig. 2. — H band intensity image of the 3x3 fields and polarization vectors (orange lines) of the polarized sources in the 
H band. The length of the vectors indicates the polarization degree. 
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Fig. 3. — (a) Histograms of polarization degree of the sources in nine fields with P > 3o\p for the J, H, and K s bands, (b) 
Histograms of polarization angle of the sources in nine fields with P > 3ap for the J, H , and K s bands. The grey dotted 
line on the histogram of the all sources show the results of Gaussian fits. The central angle and the sigma of the Gaussian 
function are listed in Table 2. 




Fig. 4. — Correlation graph of matched data in the H band. Polarization degree (left) and angle (right) are shown. The 
solid lines are linear fits with slopes of 1.52 (left) and 1.05 (right). 
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Fig. 5. — Color-magnitude (J — _ftT s vs. AT S ) diagram of our sources. Each group is indicated with a different color box: 
Group A (orange box), Group EFG (blue box), Group B (red box), Group I (green box), and Group JK (violet box). Our 
photometry data can detect more faint sources (m,K s < 16 mag) than the previous data of Nikolaev & Weinberg (2000) and 
most of the sources are distributed in the range of 14 < m,K s < 16 mag. The red and blue dashed boxes are extrapolated 
regions from Group B and Group EFG. 




Fig. 6. — (a) Histograms of the polarization degree of the sources which are classified to the Group B for the J, H, and 
K 3 bands, (b) Histograms of the polarization degree of the sources which are classified to the Group EFG and P > 3o~p 
for the J, H, and K a bands. 



144 



J. KIM ET AL. 



100 



80 



60 



40 



20 







. (o) Group B 
































. n Fft: 










i 




"-"4'"'i""T" ! 

■ - — X LT 













50 100 150 




Ol , , , , I , , , , I , , , , l_ 

50 100 150 

8 (d»g) 



Fig. 7. — (a) Histograms of the polarization angle of the sources which belong to the Group B for the J, H, and K a bands, 
(b) Histograms of the polarization angle of the sources which belong to the Group EFG and P > 3<rp for the J, H, and 
K B bands. 




SPMbi) 

Fig. 8. — Plots of the proper motion and of its error. The blue region (S/N > 2) may contain most of the foreground 
sources, and the main LMC sources are mostly distributed in the red region (S/N < 2). From the S/N ratio of the proper 
motion we try to separate sources. 
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Fig. 9. — (a) Histograms of the polarization degree of the sources which have proper motion /i > 2cr M and P > 3erp for the 
J, H , and K 3 bands, (b) Histograms of the polarization degree of the sources which have proper motion fi < 2a l _ l and P > 
3a p for the J, H , and K 3 bands. 




Fig. 10. — (a) Histograms of the polarization angle of the sources which have proper motion fi > 2cr M and P > Sap for the 
J, H, and K s bands, (b) Histograms of the polarization angle of the sources which have proper motion fi < 2a ^ and P > 
3ap for the J, H, and K 3 bands. 



